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ABSTRACT 

Dust-Obscured galaxies (DOGs) are bright 24/tm-selected sources with extreme obscu¬ 
ration at optical wavelengths. They are typically characterized by a rising power-law 
continuum of hot dust (T^i ^ 200-1000 K) in the near-IR indicating that their mid-IR 
luminosity is dominated by an an active galactic nucleus (AON). DOGs with a fainter 
24 /im flux display a stellar bump in the near-IR and their mid-IR luminosity appears 
to be mainly powered by dusty star formation. Alternatively, it may be that the mid- 
IR emission arising from AGN activity is dominant but the torus is sufficiently opaque 
to make the near-IR emission from the AGN negligible with respect to the emission 
from the host component. In an effort to characterize the astrophysical nature of the 
processes responsible for the IR emission in DOGs, this paper exploits Herschel data 
(PAGS -I- SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive 
a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust 
temperatures and masses) based on SED fitting. Of particular interest are the 24 p,m- 
bright DOGs (P 2 ig.m >lmJy). They present bluer far-IR/mid-IR colors than the rest 
of the sample, unveiling the potential presence of an AGN. The AGN contribution 
to the total 8-1000/rm flux increases as a function of the rest-frame 8/im-luminosity 
irrespective of the redshift. This confirms that faint DOGs (Lg^m < 10^^ L©) are 
dominated by star-formation while brighter DOGs show a larger contribution from an 
AGN. 
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1 INTRODUCTION 

The unprecedented sensitivity and angular resolution of the 
Spitzer Space Telescope at infrared (IR) wavelengths led 
to the discovery of a new type of galaxy that is extremely 
faint in the optical (~ 22<R<27), yet bright at mid-infrared 
wavelengths ( Houck et al.|200^ Dey et al.|2008 Fiore et al. 
2008). These sources, known to as “Dust Obscured Galax¬ 
ies” (hereafter, DOGs) in reference to the cause of their 
faintness at optical wavelengths, have extremely red optical- 
to-IR colors (/^(24/rm)//^(R) > 982). The incidence of 
DOGs is relatively low: only 8% of 24 fim detected sources 
are DOGs, while ~40 % of the sources in the 2 deg^ COS¬ 
MOS held optical catalog have similar R-band magnitude 
[22.4-26.4]. However, their contribution to the total IR out¬ 
put of the Universe at 2 : ~ 2 where their source numbers 
peak is estimated to be at least 30% (Riguccini et al. ]MT| . 
This contribution increases to 50 % when considering the 
high-luminosity tail of their distribution at these redshifts 
(i.e., I/iR > 10^^ Lq ; e.g., Riguccini et al. 2011). DOGs 


and ULIRG 


class of galaxies (e.g., Dey et al. 

to 

o 

o 

00 

Buss- 

mann et al. 

2009 Riguccini et al.||2011|. Such luminosities 


require signihcant amounts of dust-heating, most probably 
arising from star-formation and/or high levels of nuclear ac¬ 
tivity (i.e., active galactic nucleus or AGN). A number of 
recent studies have split the DOG population along these 
lines: i.e., DOGs showing a “bump” at 1.6 fim indicative of 
star-formation (|Farrah et al. 2008 Desai et al. 2009 here¬ 


after bump DOGs) and DOGs displaying a rising power-law 
SED in the near- to mid-IR bands, suggesting a dominant 
AGN ( Houck et al.|[2005 Weedman et al. 2006 hereafter 


PL-DOGs). Estimating the star-formation rate of the latter 
has proved extremely difficult due to the dominant AGN 
component washing out any host galaxy signatures. 

The faintness of DOGs at optical wavelengths has made 
the characterization of their physical properties particularly 
challenging. The launch of the Herschel Space Telescope 
in 2009 provided a new window onto these galaxies that 
is largely independent of dust obscuration, thereby giving 
us the clearest view yet of these galaxies. The wavelengths 
probed by Herschel cover the peak of the spectral energy dis¬ 
tribution (hereafter, SED) of DOGs at the redshifts where 
their numbers are highest (i.e., 1.5 <2^3). This allows us 
to accurately constrain important properties, including the 
total IR luminosity as well as dust temperature and mass. 
The aim of this work is to use the combined diagnostic pow¬ 
ers of both Spitzer and Herschel observations to determine 
how these properties relate to the dominant source of energy 
in these galaxies be it AGN, intense star formation or a com¬ 
bination of both. Eor this we use a sample of Spitzer/MIPS 
24 /rm-selected DOGs (satisfying ^ 24 ^^ > 0.08 mjy) se¬ 


lected from the COSMOS held (Scoville et al.|2007l, and 


detected in all 5 Herschel bands. We calculate the contri¬ 
bution from AGN and/or star-formation to the total energy 
output of these galaxies via SED htting and relate this to 
their dust temperature and masses. 

The paper is organized as follows. Our data are de¬ 


scribed in Sect.[^ the far- to mid-IR colors of DOGs sources 
are detailed in Sect.|^ The SED-htting procedure used and 
the results obtained are described in Sect. [4] In Sect. we 
present the model and results on the dust temperature and 
mass of our DOG sample and discuss if the presence of AGN 
signatures induce a particular trend in the Tdust distribu¬ 
tion. We discuss our results and present our conclusions in 
Sect.[^ Throughout this paper we assume a ACDM cosmol¬ 
ogy with Ho=70 km s“^, Tim = 0.3, and Ga = 0.7. Unless 
otherwise specihed, magnitudes are given in the AB system. 


2 DATA 

The sample of DOG sources is selected from the deep 
Spitzer/MIPS observations of the 2 deg^ COSMOS held 
([Sanders et al.||2007|). Our starting point are the 24 prm de¬ 
tected sources from the catalogue described in |Le Flocdi] 
et al. (2009) (see also Riguccini et al.||2011 |. We note that 
other studies further require a source to satisfy / 24 ^m > 
300 mJy in order to classify it as a DOG (e.g., Dey 


et al.|[2008|. In this study we consider all sources satisfying 


fv{24:pLm)/ft,{R) > 982 as DOGs. Furthermore, DOG stud¬ 


ies focussing on heavily-obscured AGNs (e.g., Fiore et al. 
2008| 2009| Treister et al. 20091 also impose an additional 
R-K>4.5 (vega) cut. 


2.1 COSMOS observations 

COSMOS is a wide-area equatorial held with deep cover¬ 


age at all wavelengths spanning radio to X-rays (jHasinger 
et al.|2007 Schinnerer et al.|2007 Elvis et al.|2009| . Crucial 

for this study is the deep IR coverage of this held, partic¬ 
ularly at mid- to far-IR wavelengths by the Spitzer Space 
Telescope with the MIPS instrument ( Le Floc’h et al.|2009 ) 
and, more recently, with PACS ( jPoglitsch et al.||2010p and 
SPIRE (Griffin et al.|2010l onboard Hershel iPilbratt et al. 

Moj). 


The extensive UV to near-IR coverage of COSMOS 
(e.g., Taniguchi et al. 2007| [Uapak et al.j 20071 allows for 
precise photometric redshifts (hereafter, photo-z) to be de¬ 
rived for extragalactic sources within this held. For the 
photo-zs used in this work, we use an updatec j version of 


the photometric redshift catalog of Ilbert et al. ( [2009 1 that 
provides photo-zs for 1400 237 i’^-detected sources among 
the 2 017 800 sources of the COSMOS photometric catalog. 
These redshifts are obtained with an unprecedented accu¬ 
racy, with a dispersion of o'az/(i-i-z) = 0.012 for sources sat¬ 
isfying i\g <24 and z<1.25. More relevant to this study — 
where we focus on dusty 24/im-selected sources that are very 
faint at optical wavelengths — is their comparison with the 
optically-faint spectroscopic sample from the z-COSMOS 


survey (Lilly et al. 20071 where Ilbert et al. (20091 re¬ 


port a dispersion of only crA 2 /(i+ 2 ) =0.06 for sources with 
23< i'^g <25 at 1.5< 2 <3. Given their accuracy for faint 
sources, we use these photo-zs for our 24 ^m sources, match¬ 
ing their optical counterparts following the procedure out¬ 


lined in Le Floc’h et al. (2009) and Riguccini et al. (2011). 


^ Luminous Infra-red Galaxies with 10^^ Lq <L/r < 10^^ Lq 
and Ultra-Luminous Infra-red Galaxies with L/H > 10^2 Lq (e g ^ 
Sanders et al.|1988a|b I 


^ version 1.8: the main improvements compared to [ilbert et al.| 
< |2009[ | reside in relying on the median of the PDF to define the 
“best” photo-z, instead of the minimum 
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We briefly describe this approach in the following subsec¬ 
tion. 


2.2 The far-IR counterparts of 24 pm-selected 
sources 


Our 24 pm parent sample (from Le Floc’h et al. (20091 and 


Riguccini et al. (20111) contains 29 395 sources detected at 


24 pm with F 24 f_im > 80 pJy over a total area of 1.68 deg , 
which excludes regions contaminated by bright, saturated 
objects. In the interest of focussing on the sources’ star- 
formation histories, we exclude X-ray detected AGNs down 
to a flux limit of So. 5-2 feeu = 5 x 10“^® ergcrn”^ s“ 


et al. (20091. 


on AGN catalogs from Brusa et al. (2007 20101 and Salvato 


We limit our counterpart identification to 24 pm sources 
with a 3 -(t PACS detection at 100 pm and 160 pm. SPIRE 
fluxes will be used for a subset of our sample. The GOSMOS 
field was observed as part of the PAGS Evolutionary Probe 


(PEP, Lutz et al. 20111 and the Herschel Multi-tiered Ex- 
tragalactic Survey (HerMES, Oliver et al.||2012 ) campaigns 
(i.e., PACS 100 & 160 pm and SPIRE 250, 350 & 500 pm re¬ 
spectively). The catalogs provided by PEP and HerMES cal¬ 
culate source fluxes in each of these 5 bands by performing 
PSF fitting at the positions of the 24 pm-detected sources 


from Le Floc’h et al. (20091. One of the key benefits of us¬ 


ing such 24 pm “priors” as opposed to generating blind cat¬ 
alogues, is that it helps with deblending, which is particu¬ 
larly problematic at the longer Herschel wavelengths. The 
HerMES catalog was built following the method presented 
in |Roseboom et M](|2010 |, based on the 24 pm position pri¬ 
ors from Le Floc’h et al. (20091. The PEP catalog was ob¬ 
tained using the same 24 pm priors (Berta et al. pm| . The 
reliability and the completeness of the PACS and SPIRE 


COSMOS catalogs are detailed in Lutz et al. (20111 and 


Oliver et al. (2012 I, respectively. We identify a total of 6 029 


24 pm-detected sources with a >3 -(t detection in the PACS- 
bands with Fioo^m > 3.1 mjy and FieoM^ > 6.3 mjy (see 
Tabl e[l|. We mat ch these to the catalogue of optical sources 
from |Ilbert et al.] ( |2009[ ) in order to obtain their photometric 
redshifts. 

Given the much higher density of sources detected at 
optical wavelengths in COSMOS ( |Capak et al.|[2007| ) com¬ 
pared to those detected with MIPS, a direct cross-correlation 
between the 24pm-selected catalog and the optical observa¬ 
tions could lead to a large number of spurious associations 
with optically-detected galaxies randomly-aligned close to 
the line of sight of the MIPS sources. To minimize this, 
we first matched our 24pm catalogue to the K-band cat¬ 


alogue of McCracken et al. (20101, employing a matching 


radius of 2” and following the same procedure described in 
|Le Floc’h et al.| ( |2009[ ) and |Riguccini et al.| ( |2011[ ). Of the 
6 029 sources in our sample, 5858 were found to have a K- 
band counterpart. In an attempt to reduce the number of 
non-matches, we also matched our 24pm catalogue to the 


IRAC-3.6 pm catalog from Capak et al. (20071, adopting the 


same 2” matching radius. This led to 34 additional matches, 
increasing to 5892 the number of MIPS-24 pm-|-Herschel 
sources with either a K-band or an IRAC-3.6 pm counter¬ 
part. These 5892 sources were then matched to the updated 
version of the i"'"-band selected catalog of photometric red- 


shifts from Ilbert et al. (20091 using a matching radius of 



photo —z 


Figure 1. Photometric redshift distribution of the 24pm sources 
from |Le Floc’h et al.| l [^09^ in black and of the PACS-DOGs from 
this work in purple. 


1”. Of the 5892 sources, 5768 had i’^-band counterparts and 
associated photometric redshifts, leaving 261 sources among 
the 6 029 24 pm-|-Herschel sources (i.e., « 4%) without pho¬ 
tometric redshifts. These sources were excluded from any 
further analyses. 


2.3 The PACS-DOGs sample 


Our DOGs sample was selected from the 5 892 sources se¬ 
lected at 24 pm with a 3 -(t detection in at least one of 
the PACS bands (100 or 160 pm). The DOGs criterion in¬ 
troduced by Dey et al. (20081 is based on the following: 


F’ 24 /jm/F’fl >982 and ^’ 24 ^™ >300 pJy, where the latter is 
a direct consequence of the depth of the MIPS imaging in 
the Bootes field. Considering that the source-extraction per¬ 
formed by |Le Floc’h et al.| ( |2009[ ) reaches a completeness of 
~90% with F 24 /jm >80 pJy, we extend the DOGs 24 pm- 
flux cut down to 80 pJy. The R-band magnitudes used in 
this work are from |Ilbert et al.] (|2009[) based on observations 
with the Subaru telescope by Capak et al. (20071; these in¬ 


clude a correction for Galactic extinction - not applied in 
|Capak et al.| (|2007|) - and reach a limiting magnitude of 
Mh >17.^ 

The sample contains 57 and 138 DOGs detected in 
only one of the PACS bands at 100 pm and 160 pm, re¬ 
spectively, while 119 DOGs are detected in both PACS 
bands. This amounts to a total of 314 PACS-detected DOG 
sources (cf Table[^ with M_r > 23.4, Fioo/jm >3.4 mJy and 
Fieo/jm >7.8 mJy within the redshift range 1.5< z <3 (see 
Fig- a , where the DOGs criterion is the most efficient (e.g.. 


Dey et al.|2008|[Bussmann et al.|2009[ [^12[ [Riguccini et al.| 


20111. We will use this sample in Sec.|^ 
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Figure 2. Distribution of 24 fim flux for the DOG parent sam¬ 
ple ('^2100 sources) from [Riguccini et alT] ( |2011| | (green), PL- 
DOGs (blue) and the Herschel-DOGs from this work (red). 
The Herschel-DOGs distribution peaks at higher 24 pm fluxes 
(~0.36mJy according to a Gaussian fit), compared to that of the 
whole DOGs sample (~0.14mJy). The distribution of the PL- 
DOGs selected from [Riguccini et al?] l [2011[ | also peaks slightly 
higher (~0.22mJy) than the whole DOGs population distribu¬ 
tion and is also more inclined to select 24 pm bright sources. See 
Sect.lOl for details. 



Figure 3. Distribution of 24 flux of DOGs detected in at 
least one PACS-band (314 sources) in blue and the same Herschel- 
DOGs distribution than in Fig.in red. 


2.4 The Herschelr-'DOGs sample 

Considering that we seek to undertake SED-fitting across 
the mid- to far-IR wavelength range, we further dehne a 
sub-sample of PACS-DOGs detected in all 5 Herschel bands 
; this allows for a better constraint on the peak of the DOGs’ 
SEDs. To achieve this we matched the PACS-DOGs sample 


with the SPIRE catalog ( |Roseboorn_ et ah 20101. This re¬ 
sults in 95 Herschel-detected dogR i.e. detected in the 5 
Herschel bands (see Table[^ . 

The fluxes in the PACS bands are obtained for all 
sources with a > 3-a detection. Although the SPIRE catalog 
reaches a S-cr limit of ~ 10 mjy, ~ 12 mjy and ~ 15 mjy at 
250, 350 and 500/rm respectively, the corresponding 3-a ex- 
tragalactic confusion limits are 14.4 mJy, 16.5 mJy and 18.3 
mJy ( Nguyen et al.||2010 |. In the SED-fitting procedure we 
are cautious ( [Magnelli et al.||2012a[ ) when including fluxes 
that are lower than the 3-a extragalactic confusion limits 
and use them merely as upper limits. 

Among our sample of 95 Herschel-detected DOGs, 40 
have their fluxes above the 3-a threshold only for the 250 and 
350 bands, 20 merely for the 250/im band, one DOG for 
the 350 and 500 fj,m bands and another DOG solely for the 
350 fim; 9 of the Herschel-detected DOGs have all SPIRE 
fluxes below the S-cr threshold. We quote upper limits for 
all of these cases. Only 24 sources have fluxes above the 3-cr 
limit in the three SPIRE bands. 

We acknowledge that imposing a detection in the 5 Her¬ 
schel bands will impart a bias towards the brightest and red¬ 
dest IR sources in our sample. This is shown in Fig.[^ where 
the distribution of Herschel-DOGs peaks at higher 24 /rm 
fluxes than that of the DOGs parent sample from|Riguccini| 


et al. (2011 Of particular interest is to note that although 

<0.4 mJy) are missed by the 
>0.4 mJy the Herschel- 


the faintest DOGs (F24 m>i 
Herschel-selection, beyond F 24 ^ri 
DOGs distribution is very similar to that of the DOG par¬ 
ent population (see Fig.j^. On the other hand, PL-DOGs - 
known to be mainly AGN dominated (e.g.,|Bussmann et al.| 


20091 - present a significantly stronger bias towards 24^m- 


bright sources: not only they have a 24pm flux distribution 
that peaks slightly higher than the whole DOGs population 
distribution, but their selection represents 60% of the DOG 
population with F 24 ,jm >lmjy, compared to merely a 10% 
at F 24 fj,m =0.3mJy |Bussmann et al!| ( |2009[ ). 

In this paper, two samples of DOGs are used. To study 
the IR colors of DOGs (Section|^ we use only PACS data 
and thus base our analysis on the 314 PACS-DOGs, in an ef¬ 
fort to improve our statistics. For the remainder of our study 
we restrict our analysis to the 95 Herschel-detected DOGs, 
noting that both samples probe the same DOG population, 
as illustrated by their 24 pm flux distributions on Fig.|^ We 
focus our study on the differences observed between mid-IR 
bright DOGs (F 24 ^m > ImJy) and DOGs with more mod¬ 
erate fluxes (~ 0.2mJy<F24Mm < ImJy). 


2.5 The reliability of the photometric redshifts 
for the DOGs sample 

The high accuracy of the photometric redshifts for sources 
of the GOSMOS catalogue (see Sect. |2.l| ) make them highly 
reliable for statistical studies on large (i.e., >2 000) sam¬ 
ples (e.g., Riguccini et al. 20111. However, because in this 
study we focus on significantly smaller numbers, we require 
particularly robust redshift measurements for each source. 
To ensure this, we checked the distribution of the probabil¬ 
ity density function (PDF) of the photometric redshift for 


® We note that 24 pm-selected sources with no R-band detection 

are also considered DOG sources. ^ All 95 Herschel-detected have R-band detections 
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Table 1. Number of sources 


6 029 

24 /im-sources with a 3-a PACS detection (at 100 fim 
and 160 fim) 

5 892 

sources from the previous sample with a photo-z 

314 

DOGs (i.e., F24,ij,m./Fpi >982) with ^ 24 ^™ >80/rJy 
and 1.5< 2 : <3 and with a 3-cr detection in one or 
the two PACS-bands 

Sample used for the far-IR/mid-IR color analysis 
[Sect.j^ 

95 

DOGs with F 24 ^m >80/iJy and 1.5< 2 <3 and with a 
3-cr detection in the 2 PACS bands and with a detec¬ 
tion (potentially > 3 cr) in the 3 SPIRE bands 

Sample used for the remainder of the paper 


As an initial, crude assessment of the dominant pro¬ 
cess responsible for producing most of the IR output in 
DOGs (i.e., AGN vs. star formation), we first consider the 
far- to mid-IR colors (hereafter FIR/MIR) of our sample 
(e.g., Mullaney et al.|2012 K Fig •l^shows the 100 /rm/24 fim 
and 160 /rm/24 ^m color distributions for our sample of 
PACS-detected DOGs as a function of redshift; we include 
all 24 fim -detected GOSMOS sources for comparison. We 
see no noticeable trend for the DOGs sample at z^2; the 
curves shown for the 100/24 and 160/24 median color evo¬ 
lution with redshift seem to follow the same evolution than 
that of the whole 24 ^m-detected sample. However, at lower 
redshifts the DOGs display a steeper evolution, with bluer 
100/24 colors than the non-DOG 24pm-detected sources. 


each source and divided our sample into three categories, 
according to the photometric-redshift reliability. The cat¬ 
egories are the following: (1) sources have a single, secure 
photo-z, for which the PDF has a gaussian shape with a 
single peak (39 sources); (2) those with multiple potential 
photo-zs, for which either the PDF’s peak is spread over a 
wider range of redshifts {Az ~ 0.2) or it includes a lesser 
peak which may correspond to another photometric red¬ 
shift (24 sources); (3) those flagged as presenting inaccurate 
photo-zs, because their PDF shows clear multiple peaks of 
similar strength (32 sources). For sources in the last two cat¬ 
egories, the photo-z is set to the highest peak value of the 
z-distribution and in the case of multiple peaks we keep the 
value of the subsequent peaks as secondary options. To get 
the most reliable and accurate fit to the SEDs of our sources, 
we use all of these potential photo-z values. We make a spe¬ 
cial note that six DOGs in our sample have a confirmed 
COSMOS spectroscopic redshift as part of the Fiber Multi 
Object Spectrograph (FMOS) spectroscopic redshift catalog 
(Kartaltepe et ah, in prep.) which we use in our analysis for 
higher accuracy. 


3 RESULTS: FAR-IR/MID-IR COLORS OF 
EXTREMELY MID-IR BRIGHT DOG 
SOURCES 


The DOG sources are not only an extreme sub-sample of 
ULIRGs but also represent a mix between sources domi¬ 
nated by star-formation and those dominated by AGN ac¬ 
tivity (e.g., Houck et al.|2005 Fiore et al.|2008 2009 Buss- 


mann et al.||2009 Melbourne et al.||2012 |. In this paper we 

seek to quantify the AGN contribution of these sources and 
study the evolution of this contribution with respect to other 
galaxy properties, including redshift, the 8^m rest-frame 
luminosity, total IR luminosity, dust temperature and dust 
mass. 

Studies in the past years have explored the PL- and 
bump-DOGs population (e.g., Pope et al.||2008 Melbourne 
et al. 20091. It has been well established that PL-DOGs have 


an AGN contribution to their near-IR emission and that 
their far-IR emission is most likely dominated by star for¬ 
mation) CaIano££r^aIQ20^ I. In this paper we aim to gauge 
the AGN contribution of these DOGs using Herschel far-IR 
data. 


We find that the FIR/MIR distribution of both the bulk 
of the 24 fim comparison sample as well as the majority of 
our DOGs are well represented by the star forming tem¬ 
plates from (Ghary & Elbaz|2001| hereafter GEOl) with IR 
luminosities L7R^=l(r^^^^°L0. This is to be expected 
given that IR-selected galaxies at z>l tend to be of the 
LIRG or U LIRG class (e.g, |Le Floc’h etH][200^ [Magnelli 


et al. 20091. We emphasize that although recent work has 


shown that GEOl local ULIRG SEDs are not good fits to 
z~2 star-forming galaxies with similar IR luminosities (e.g., 
Elbaz et al.(Mol[mT|[Nordon et aI.|2010|[M^, the GEOl 


ULIRGs templates are good fits to our DOGs; 17 of our 
sources are fit using these templates with a <5 (Sec.|^. 


We consider the particular case of the brightest DOGs 
in our sample (i.e., F 24 ^m >1 mjy) and find that they show 
a particular behavior in their FIR/MIR colors as a function 
of redshift: the brightest DOGs in our sample show signif¬ 
icantly bluer PACS/24 ^m colors than the general 24 ^m- 
detected population (i.e, with F24Mm >80 uD)- We discard 
the possibility that a variation in the Photodissociation re¬ 
gions (PDR) component and/or variation in the intensity 
of the field is responsible for the bluer color of these bright 
DOGs, by comparing to the 100 ^m/24 fim colors derived 
from the templates of |Magdis et al.| ( |2012[ ) (see Fig.[^. These 
SED templates are based on stacked ensembles at different 
redshift intervals, considering the varying radiation field and 
PDR contribution to ULIRGs as a function of redshift; for 
our study we rely on their starburst-dominated templates 
at the two relevant redshift intervals: 1.75 < z < 2.25 and 
2.27 < z < 3.0. We compare to the FIR/MIR colors of 
AGN/galaxy composites - using the intrinsic AGN SED of 


Mullaney et al. (20111 and assuming different AGN contri¬ 


butions (25, 50, 100%) at 100 /rm and 160 /rm - and find 
significant similarities, suggesting that the brightest DOGs 
have a significant AGN contribution. As we consider higher 
redshifts, the median FIR/MIR color of these bright DOGs 
point towards a lower fraction of the AGN contribution, con¬ 
sistent with the SFRs of galaxies of similar mass increasing 
with redshift (e.g., [Brinchmann et ar]|2004[ |Daddi et al.| 


2007 Pannella et al. 2009 [Magdis et al. 2010[ |. We check 
the validity of these trends in the following section by look¬ 
ing at the AGN contribution (based on SED-fitting) as a 
function of the 24 /rm flux. 


© 2002 RAS, MNRAS 000,[iHT7l 



































































6 Riguccini et al. 




Figure 4. Left: Distribution of lOO/rm/24/rm color as a function of redshift for all the 24 /rm-selected sources detected at 100 fim (black 
open circles), the 176 DOGs detected at 100 fim with -F 24 (im >80 /rmjy (grey open squares) and with F 2 A)im >1 rnjy (green filled circles). 
The black, grey and green solid curves, from top to bottom, represent the median of all the 24 /rm-selected sources with F 24 ,fxm >80 
pjy, the median of the 24 pm-sources selected as DOGs with F 2 Aiim >80 /rmjy and the median of the brightest DOGs sources (i.e., 
with F 24 ^m >1 uijy), respectively. Errors on the median are calculated as quadratic propagation of uncertainty. We show the expected 
flux ratios (red dashed-triple dotted tracks) in the case that 25, 50 and 100% of the flux in the 100 /rm band result from an AGN 
component l |Mullaney et al.|2011[[2012[ |. We also include for comparison the star-forming ULIRG GEOl templates for an IR luminosity 
of 10^^ Lq (bottom blue-dashed line) and 10^^'® Lq (top blue-dashed line), as well as the template from [Magdis et al.| ( |20i^ (light 
blue dotted-dashed line). The observed PAGS/24 colors of the bulk of the 24 /im sources and that of the DOG sources are consistent 
with these ULIRG templates. Right: Distribution of 160/4m/24/rm color as a function of the redshift for all the 24 /rm-selected sources 
detected at 160 fira (black open circles) and the 257 DOGs detected at 160 fim. The colors and tracks are the same as on the left panel. 
For clarity, we do not overplot the GEOl templates on the right panel as they would give the same results than for the 100/24 color. 


4 AGN CONTRIBUTION TO THE TOTAL 
REST-FRAME S-lOOO/iM LUMINOSITY IN 
DOGS 

Based on FIR/MIR colors, bright DOGs likely contain an 
AGN component, contributing partly or even dominating 
their IR luminosity. To have a better understanding of these 
sources and to have a global view of their stage in the evolu¬ 
tion of the galaxies, it becomes important to know the exact 
contribution of a potential AGN to their total S-lOOO/rm rest 
frame luminosity. In this section we present our method to 
determine the potential contribution of an AGN component 
to these DOG sources and show our results on the variation 
in AGN contribution with the 24 /rm flux. 


4.1 Method: SED-fitting procedure 

Studying the SED of a galaxy provides insights to the physi¬ 
cal nature of the underlying continuum source and can unveil 
the presence of an AGN. The impact that an AGN contribu¬ 
tion has on the shape of the SED is distinct from that of dust 
heated by star-forming activity. However, deriving the SED 
of a galaxy is not an easy exercise especially in the case of an 
AGN where the imprints of the host galaxy is always present. 
In our study we use the IDL-based SED-fitting procedure 
DecompIR, detailed in [Mullaney et ST] ( |20li| . Combining 
a set of five starburst templates and an average AGN tem¬ 
plate, this approach is aimed at fitting the IR photometry of 
composite galaxies and to measure the AGN contribution to 
their total IR output. A method is used to know which 


combination of these templates best fits the data; i.e., the 
combination with the lowest associated value is adopted 
as the best fit. The validity of this procedure as an accurate 
way to determine the AGN contribution to the total IR out¬ 
put of composite galaxies has been verified by several tests 
lead by Mullaney et al. (20111, including a comparison with 
alternative measures of the AGN contribution (e.g., emission 
line diagnostics). Although there are significant uncertain¬ 
ties associated to the precise AGN contribution to an indi¬ 
vidual galaxy, this approach is adequate from a statistical 
point of view (i.e., large samples, average SEDs). 


We apply the DecompIR procedure to our sample 
of DOGs sources. Nonetheless, considering that DOGs 
have ULIRG-class luminosities (e.g., Bussmann et al.||2009 


Riguccini et al.[2011|, we add two ULIRGs templates with 


Liu ~ 10 Lq and Lm = 10 “Lq to the set of starburst 
templates from Mullaney et al. (20111 to fully cover the lu¬ 
minosity range of our sample. The ULIRG templates are 
taken from GEOl as they build their library from SEDs and 
models that only take into account star formation activity. 
The median PACS/MIPS colors of the bulk of the DOGs 
sample are well represented by these ULIRGs templates (see 
Fig.|4f, motivating their use as part of our SED-fitting pro¬ 
cedure. To determine the AGN contribution to the IR lu¬ 
minosity of our sources, we first derive a best SED-fit with 
templates based only on star-forming sources. If the star¬ 
forming templates do not provide a satisfactory result, an 
AGN component is added and the SED-fitting proceeds with 
a composite spectra. We consider the AGN component as a 
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reasonable option only if it improves the of tho fit by at 
least 50%. 

We implement our SED-fitting procedure to each DOG 
source in our sample. Each of the 95 sources are first fit 
by a star-forming component only and then by a composite 
spectra when the x^ from the star-forming fit is > 20. In 
the case of sources with a less-accurate photo-z, the source 
is fit with a star-forming template for all possible photo- 
zs obtained from the PDF (see Sect. 2.51. If none of these 
fits are suitable, we add an AGN component and the fits for 
each possible photo-z are performed once again. Our method 
is robust in fitting most of our sources (90%). Out of 95 
DOGs, 71 are fit with a star-forming galaxy template only, 
15 require an AGN component to the fit, and in 9 cases, 
no reliable fit was obtained either using a star-forming-only 
template nor a composite spectra. The failure of successfully 
fitting these sources could be due to the fact that we cannot 
reproduce the SED of these sources; this is most likely due to 
a wrong redshift, even after probing the different possibilities 
indicated by the PDF. 

Due to the expected uncertainties on the AGN frac¬ 
tion for an individual galaxy, we implement the SED-fitting 
procedure to average DOG SEDs. We divide our DOG sam¬ 
ple in 3 different redshift bins (1.5<z<2.0, 2.0<z<2.5 and 
2.5<z<3.0) and in four 24 /rm-fiux bins (0.09< F 24 ,jm <0.24, 
0.24< F24 m™ <0.65, 0.65< F24/jm <1.76, and 1.76< 
F 24 /jm <4.74mJy). For each redshift bin and each 24 /rm- 
flux bin, an averaged SED is calculated at 8 fim, 24 fim, 
100 fim, 160 /rm, 250 350 and 500 /rm, leading to a to¬ 

tal of 12 average SEDs. The results of the fits are shown in 
Fig[^for the average SEDs and in Fig[^&[^for the individual 
sources fitted with an AGN. 

The results from the SED-fitting procedure on the av¬ 
erage DOG SEDs are presented on Table[^ including the 
AGN contribution to the total IR flux (fAGiv), the best- 
fit template and corresponding x^- Our method to get the 
AGN contribution does not appear to be biased toward one 
specific template. Independently of the AGN contribution, 
the two templates that were the most successful at fitting 
the average SEDs were the CEOl templates for IR lumi¬ 
nosities of 10^^ Lq (CE12) and 10^^'® Lq (GE12.5). This is 
as expected, since these templates have PAGS/MIR colors 
consistent with that our DOG sources. 


4.2 Results 

We list the mid-to-far IR photometry and the total IR lumi¬ 
nosity of the 95 Herschel-DOGs on Table 3, indicating also 
whether an AGN component is included as part of the SED 
fit. The redshift distribution of our DOG sample peaks at 
z~2, allowing us to use the 24 /rm band as a probe of the 
mid-IR emission close to the rest frame 8 /rm. This allows 
us to derive a Lg/jm that minimizes the dependence on the 
choice of SED template used to perform the fc-corrections. 
To determine the rest-frame Lg^m we interpolate the CEOl 
library at the redshift and flux of each average SED. We are 
then in a position to study the fractional AGN contribution 
to the total 8-1000 /rm output as a function of the 8fj,m rest- 
frame luminosity (Lg^m); we do this for the average SEDs at 
the three redshift bins: z = 1.5 — 2.0, 2.0 — 2.5 and 2.5 — 3.0 
(see Fig. 1^. 

We find that the AGN contribution increases globally 


with increasing Lg^m for all our redshift bins (see Fig. [^. 
This confirms the findings of Pope et al. (20081, where they 
report - based on mid-IR colors of 79 sources within the 
GOODS field and with 24 /xm fluxes down to 100 /rjy - 
that low-luminosity DOGs are primarily powered by star- 
formation activity. However, only MIPS-70/im observations 
were available for their analysis and the inability to sam¬ 
ple properly the peak of the SED lead to large uncertainties 
on the derivation of the dust temperature and the AGN 
contribution; Penner et al. (20121 extend the study out to 
far-IR wavebands but miss the faintest DOGs by focussing 
on GOODS DOGs with luminosities 10^^ Lq < L/h < 10^® 
Lq. Furthermore, Fiore et al. (20081 claimed that even faint 
DOGs show evidence of hard X-ray emission, suggesting the 
presence of an underlying AGN contribution. Within this 
context our work - based on the large 2 deg^ area of the 
COSMOS field and the good sensitivity of the MIPS-24 fim 
observations that allows us to sample a wider range of 24 fj.m 
fluxes (i.e, 80 /rJy < F 24 ,iim < 5 mjy) and the access of far-IR 
data with Herschel - allows us to conclude that faint DOGs 
are mainly star-forming systems while brighter sources be¬ 
come dominated by an AGN. 

By separating our sample in redshift bins. Fig. also 
shows that the relation between AGN fraction and Lg^m 
evolves with redshift: the slope of the AGN contribution 
with respect to the Lg^m, is steeper at low redshifts, while 
at higher redshifts the AGN contribution is less important. 
This is consistent with the results from IMerloni fc Heiiizl 
(20081, where they find that although the accretion rate 
density onto supermassive black holes (SMBH) and star- 
formation rate (SFR) densities increase from 2 ~ 0 to then 
decrease beyond z~2, the decrease in SMBH activity is 
sharper than that of the SFR. We note that the uncertainties 
on the derived AGN contributions are calculated from the 
formal error output resulting from the x^ in the SED-fitting 
procedure. 

In the interest of studying the star formation activity 
in our sources, we use the results from our SED decom¬ 
position to extract the AGN contribution and calculate IR 
luminosities only due to star formation. The resulting values 
span the range of 10 ^^ Lq < Lm < 10^^ L©, correspond¬ 
ing to one order of magnitude fainter than the analysis by 


Penner et al. (20121. We study these IR luminosities as a 


function of the rest-frame 8/rm luminosity and find that for 
a given 8/rm luminosity DOG sources whose SEDs are best 
fit with the addition of an AGN component exhibit signifi¬ 
cantly lower IR luminosities than DOGs fit with a host-only 
component (see Fig.j^ top panel). In fact, for a given 8/rm 
luminosity the majority (~75%) of DOGs fit by a host-only 
component display similar IR luminosities to the median 
star-forming galaxies within the GOODS-Herschel sample 
from [Elbaz et a~ (20111, while AGN-DOGs populate the 
lower-tail of IR-to-8/im luminosity ratios (IRS = L/rj/Lg^m; 
see Fig.j^ bottom panel). We observe an anti-correlation 
between the IRS ratio and the Lg^m with the AGN-DOGs 
populating the brightest Lg^m-end. 

For each DOG source in our sample we convert the 
IR luminosity into SFR ac cording to |Kennicutt| ( |1998| ) and 
adopt the stellar mass from Ilbert et al. ( 2009[ ). We consider 
the redshift evolution of the specific SFR (sSFR = SFR/M*) 
of our DOG sample and find that the majority of DOGs 
with no AGN component display sSFRs that place them at 
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Figure 5. SEDs of the 12 average templates of DOGs galaxies from our sample. The solid line corresponds to the total SED fit, the 
dashed line is the host template and the dotted line is the AGN component. The name of the template used is written on each panels. 
The flux bins are specified on the left side and the redshift bins are written on the top. 
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Table 2. Results from the SED-fitting for the average SEDs per redshift bins and per flux bins. Jagn is the percentage of total 8- 
1000/rm flux that comes from AGN. We specify the star-forming template best fitting the host galaxy (SB for the starburst templates 
from | |Mullaney et al.|2011[ | and CE 12.5 and CE 12 for the CEOl ULIRG templates. 


Redshift 

Flux 

Jagn 

error 

Template 


Lir 


(mJy) 

% 

% 

(Host galaxy) 


L© 

1.5<z<2.0 

0.09<f<0.24 

0 

- 

CE 12.5 

4.4 

9.60e-|-ll 


0.24<f<0.65 

0 

- 

CE 12 

5.3 

1.29e-(-12 


0.65<f<1.76 

0 

- 

SB3 

5.3 

1.69e-(-12 


1.76<f<4.74 

83 

3.5 

SB2 

6.5 

2A2e+12 

2<z<2.5 

0.09<f<0.24 

0 

- 

CE 12.5 

6.9 

1.50e-(-12 


0.24<f<0.65 

0 

- 

CE 12.5 

34.6 

3.05e+l 


0.65<f<1.76 

17 

5.6 

CE 12.5 

0.8 

5.52e-(-ll 


1.76<f<4.74 

86 

15.7 

CE 12 

0.8 

2.43e-(-12 

2.5<z<3.0 

0.09<f<0.24 

0 

- 

CE 12.5 

21 

2.42e-|-12 


0.24<f<0.65 

20 

3.5 

CE 12.5 

0.4 

5.60e-|-ll 


0.65<f<1.76 

38 

3.2 

CE 12.5 

0.1 

1.69e-(-12 


1.76<f<4.74 

59 

4.5 

CE 12 

13 

6.68e-|-12 



Figure 8. Evolution of the contribution of the AGN component 
to the total rest-frame 8-1000/rm flux of the sources as a function 
of the 8/rm rest-frame luminosity (Ls^jm). The AGN fraction is 
given for 3 different redshift bins: 1.5<z<2.0 (solid black line), 
2.0<z<2.5 (pink dot-dashed line) and 2.5<z<3.0 (bleu dashed 
line). The AGN contribution is obtained from the fitting proce¬ 
dure described in Sect. on the 12 average SED (see Table[^. 
The general trend is an increasing contribution of the AGN com¬ 
ponent with respect to the Ls/jm of the source irrespective of the 
redshift range. 


or above the main sequence (MS) from Elbaz et al. (20111, 
while 50% of the AGN-DOGs show significantly lower sSFR 
values (i.e., they lie below the MS)(see Fig. |10[ |. Sources that 
lie a factor of 2 above the MS are considered as “starbursts” 
by Flbaz et al. (20111. All but three host-component galaxies 
lie within a factor of 2 around the MS or in the starburst’s 
zone. The distribution in sSFRs shown in Fig. 1 10 1 highlights 
the composite nature of the DOG population: some DOGs 
are dominated by starburst activity, the majority is undergo¬ 
ing star-formation as part of the MS, while others are domi¬ 


nated by an AGN. This prompts the idea that DOGs are at 
the crossroads of the ULIRG-quasar scenario proposed by 
( Sanders et al.|jl988a|b Bussmann et al.||2012 (, with AGN- 
DOGs being closer to a quasar phase, where the AGN has 
already started to quench the star formation (explaining the 
lower sSFR observed on Fig. 101. 


4.3 Comparison with IRAC-color AGN selection 
criteria 


Our SED-fitting analysis identifies 15 Herschel-detected 
DOGs with an important AGN contribution to the total 
IR output. We compare our AGN classification of DOGs to 
prior approaches relying on an IRAG-color selection. Fig. 


11 


shows the IRAG-color selection of AGNs by [Lacy et al. 
(20041, as well as the refined IRAG-color selection of Donley 


et al. (20121, which also includes a power-law criteria in the 


mid-IR: Ss.e < S4.5 and S4.5 < Ss.s and S5.8 < Ss.o- 

The majority of our sources display IRAG colors con¬ 


sistent with the criteria of Lacy et al. (20041, which would 


suggest that 90% of our DOGs are AGNs. However, our 
SED-fitting analysis indicates that only ~15% of our sources 
have a large AGN contribution. Based on this we conclude 


that relying on the AGN criteria of Lacy et al. (20041 would 


lead to a lack of precision in selecting AGNs versus galaxies 
dominated by star formation. On the other hand, more than 


50% of our AGN-DOGs lie within the Donley et al. (20121 


criterion, suggesting that it is a more reliable way of select¬ 
ing AGNs in DOGs when considering merely IRAC colors. 
However, from the 19 DOGs that lie within the AGN-criteria 
of [Donley et al.| ( |2012| ) - and excluding the 3 that do not 
follow the power-law criteria required by the authors - only 
9 are classified as AGNs following our SED-fitting analysis. 
That is, 40% of the Herschel-DOGs with IRAC colors con¬ 


sistent with the criterion of Donley et al. (2012 I do not have 


a significant AGN contribution according to our analysis. 
Of particular interest is that out of all our AGN-DOGs, six 
(i.e., ~40%) are not identified as AGNs based on the crite- 
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Figure 6. Following the same format as Fig.[^ this figure shows the results of the SED-fitting procedure for DOG80, one of the 15 
DOGs which require an AGN component: the first 7 panels are the results of the SED-fitting with a host component only and the last 
panel (bottom middle panel) is the acceptable fit, with a contribution of an AGN component (20% <fAGN <40%). 


ria by Donley et al. (20121, four of which do not follow the 


power-law criterion required. 


On the one hand, our AGN classification is based on 
the availability of far-IR data for obscured sources such as 
DOGs. On the other, |Lacy et al.| ( |2004[ ) and [Donley et al.| 
(20121 classify sources as AGN-dominated based on IRAC- 


color selections. When considering these selections side by 
side, we draw two main conclusions: (1) non PL-DOGs po¬ 
tentially host an AGN that may dominate the far-IR regime 
even when missed by the IRAG-color selection criteria of 
[Lacy et al.j ( |2004| ) and [Donley et al.j ( [2012[ ); and (2) PL- 
DOGs with an AGN according to our SED-fitting procedure 
can be missed by IRAG colors criteria. We conclude that 
our method provides an alternate means of determining the 
composite nature of DOGs. 


5 DUST TEMPERATURES AND MASSES 

It has been well established that interstellar dust absorbs a 
large fraction of the UV/optical radiation from DOGs and 
reemits it in the IR ( jPenner et al.|20T2 1 As such, it is essen¬ 
tial that we understand the dust properties of these galaxies 
if we are to understand this potentially-important popula¬ 
tion of galaxies. In this section, we derive the dust temper¬ 
atures and masses for our sample of DOGs. The availability 
of far-IR data from Herschel is crucial to obtain these prop¬ 
erties. We are now able to extend previous studies on DOGs 
that did not have access to such high-quality far-IR data 
(e.g., [Dey et al.[[2008[ [Bussmann et al.[[2009[ [201^ . We are 
also in a position to compare results with other recent stud¬ 
ies using (limited) Herschel data on DOGs, including that 
of SPIRE-detected sources (down to only F 24 ;jm > 0.3 mjy) 

Melbourne[ 


with spectroscopic redshifts in the Bootes field by 


et al. (20121 and the study by Calanog et al. 


(20131 on 


SPIRE-detected DOGs within GOSMOS; no PACS data 
were available for either study. 
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Figure 7. Similar to Fig.[^ SEDs of the 15 DOGs which require the contribution of an AGN component. 
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Figure 9. Upper panel: Comparison of Ljr with Lg (rest-frame 
8 fim) for DOGs from our sample (galaxies with host component 
only are marked with black open circles and AGNs are marked 
with red filled circles). For comparison, we show the median loca¬ 
tion of star-forming galaxies from [Mbaz et ar] l |2011[ | (solid line), 
with the dashed lines showing the 68% dispersion. Lower panel: 
Variation of the IR8 (=L/i?/L8) ratio with the 8/im luminosity 
for our DOG sample, following the same color code as in the 
upper panel. For comparison we also plot galaxies from |Elbaz| 
|et al.| ( [20TT] |, including local galaxies (blue crosses), star-forming 
galaxies at z>1.5 (orange triangles) and AGNs at z>1.5 (green 
asterisk). See text for details. 


5.1 The Single Temperature Model 


The availability of Herschel far-IR data allows us to con¬ 
strain the peak of the SED in the far-IR regime and to cal¬ 
culate the dust masses and temperatures of our galaxies with 
a higher accuracy than previous studies. DecompIR does not 
provide information on the dust amount of our sources; we 
ht a blackbody spectrum of temperature T using far-IR 
data (see Amblard et al.|[2014| . To summarize, we perform 
a single temperature fit (hereafter IT model) with an emis- 
sivity, /3, of 1.5 to fluxes long ward of Xrest-frame > 40 /rm. 
The luminosity is then expressed as 

L{v) <x B{v,Tei)vR (1) 

Considering Xrest-frame > 40 pLm^ we avoid emission from 
the AGN that can boost the dust temperature and bias our 
results (Netzer et al. 2007 Mullaney et al. 20111. For this 


reason we restrict ourselves to using only the PACS 100, 160 
fim and SPIRE 250, 350, 500 ftm bands for galaxies with 



Figure 10. Redshift evolution of the specific SFR (sSFR = 
SFR/M*) of DOGs; we distinguish between DOG sources whose 
SEDs are best fit with a host-only component (black open circles) 
and with the addition of an AGN component (red filled circles). 
The solid line represents the star-forming main sequence from lEl-l 
|baz et ah] | |2011[ | and the dashed lines are a factor 2 above and 
below this fit. See text for details. 



Figure 11. IRAC colors of host-component DOGs (black open 
circles) and of the AGN-DOGs (red filled circles) from this work 
; we also indicate the IRAG colors of the 9 sources for which no 
SED fit was possible (blue crosses; see Sect. 4.3 for details). The 
solid line box is the AGN selection crite rion from [D onley et al.| 
(|2012[| and the wider dashed box is from|Lacy et al.| ||2004[|! 


z<1.7 and only the SPIRE 250, 350, 500 fim bands (when 
available) for z>1.7 galaxies. We require a minimum of 3 
data points to fit the SED. At z<1.7, 15 DOGs comply with 
this requirement (from the 16 DOGS at z<1.7). At z>1.7, 
from the 24 sources that have SPIRE fluxes above the 3-cr 
limit, 22 have 3 data points. 

We enforce the dust temperature to be constrained 
between 10 and 95 K, the luminosity between 10'^° and 
10 '" La, acknowledging the high luminosities of our sample 
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Figure 12. Distribution of the dust temperature of the 24 DOGs 
with a detection in the 3 SPIRE bands with the IT model (22 
sources at z> 1.7 and 2 sources at lower redshifts). The dust 
temperature peaks at T^^ ~40 K and is comprised within the range 
24< Tdust <65 K. 


(see F ig.[^. We observed the same definition as Amblard 


et al. (20141 for the dust luminosity and for the dust mass : 


Figure 13. Evolution of the dust temperature obtained with the 
IT model as a function of the AGN fraction. The results are ob¬ 
tained for the 12 average SEDs presented in Table[^and detailed 
in Sec. |4.1| To improve clarity, we have slightly changed the x 
values for the 6 points with no AGN fraction in order to exhibit 
more clearly the error bars on the figure. The blue line is the best 
fit with a slope of 0.33di0.08 for the 6 average SEDs with an 
AGN contribution. The dotted lines represent the 1 a error. 


Ld{\) = AiiMdK{\)B{\Td) 


( 2 ) 


5.2 Discussion 


Md = Ld/ j 4nGOB{X,T)dX 


(3) 


where k is taken at 850 /rm ([ Dmme^et;^ 2000|) and 
equal to 0.077 kg“^/ni^ (Draine & Lee||1984 Hughes et al. 

We find a median dust temperature of 
Td ~(40.6±9.2)K for our sample. The dust tempera¬ 
tures of our DOGs are in overall good agreement with 
estimates from the literature for other samples of DOGs. 


Bussmann et al. (20091 predicted a high dust temperature 


for DOGs sources with >35-60 K, but this estimate 
was mainly based on observed wavelengths shortward of 
~350/im: only 4 of their 12 DOGs have 350/im fluxes, the 
rest of the sample has only upper limits). [Melbourne et ah] 
(20121 find lower dust temperatures for their Herschel- 


detected DOG sources (i.e., 20<Td <40 K). They split 
their sample into bump DOGs and PL-DOGs and find 
that the PL-DOGs are less likely to be detected at far-IR 
wavelengths using SPIRE than the bump DOGs. They also 
claim that SPIRE detections are biased towards very cold 
sources. We note that our range of temperature (~24-65 K) 
is wider than that of Melbourne et al. ( 2012|). Our dust 
temperatures are in good agreement with [Galanog et al.j 
(2013 within uncertainties) : they find Td = (37±6)K for 


detected PL-DOGs and Td = (35 ± 7) K for detected bump 
sources. As raised in the literature (e.g., [Melbourne et ah] 
2012), using only SPIRE data tends to underestimate the 


dust temperature. Therefore, we need to be cautious in our 
analysis since more than half of our dust temperatures are 
obtained using only SPIRE data. 


5.2.1 Effect of the AGN contribution on the dust 
temperatures 

With the goal of improving our accuracy in deriving dust 
temperature and mass, we use the average SEDs obtained 
per bin of 24 /rm fluxes and per bin of redshift (see Sec. 4.1 
and in Table[^ for details) instead of the SEDs of individ¬ 
ual DOG sources. The dust temperatures have been derived 
for the average SEDs following the procedure described in 
Sec. [5.1[ These are presented in Fig.[^as a function of AGN 
fractional contribution. Roughly half of the average SEDs 
have no AGN contribution and show a wide range of Td 
as seen on the left side of Fig. [13[ The two sources with no 
AGN fraction and with the highest dust temperatures show 
extremely large errors bars on the dust temperature (i.e. 
±15 K and ±16 K) while the average dust temperature er¬ 
ror for the sample is around 9 K. The rest of the average 
SEDs have an AGN contribution ranging from 20% to al¬ 
most 90% and are within the same Td range as the sources 
with no AGN component. The DOGs’ average SEDs with 
an AGN contribution display a correlation between the AGN 
contribution and the dust temperature. We perform a best- 
fit on the 12 data points, taking into account the errors 
both on the x and y axes and find a slope of 0.33±0.08 with 
a reduced x^ of 0.94. To insure the validity of our fitting 
method, we also perform a fit T = cste with 


cste = 


E T /(j X * 
E I/ct't ^ 


= 42.6 A. 


The reduced x^ for this flat fit being 2.53, this gives us a 
strong indication that our previous fit is valid. We also es¬ 
timate the Spearman’s (rho) rank correlation of Tdust and 
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the AGN percentage. The significance is low (0.04), which 
indicates a significant correlation. The correlation between 
the dust temperature and the AGN percentage for the aver¬ 
age SEDs with AGN contribution is thus real, even though 
the slope is small. However, we do not confirm the presence 
of a general trend between the AGN fraction and the dust 
temperature of the sources since half of the data points with 
no AGN activity present similar dust temperatures than the 
average SEDs with a large contribution from an AGN. 


5.2.2 Effect of the AGN eontribution on the dust masses 


In addition to dust temperatures, the IT model fitting pro¬ 
cedure also provides us with the dust masses of our sam¬ 
ple sources. We obtain a range for the entire sample of 
7 X10^ < Mdust < 10® Mq and a median dust mass of 
~ (3 ± 3) X 10® M 0 . Our results are in very good agreement 


with Bussmann et al. (2009 \ who found a median dust mass 
of 3x10®Mq for their sample of 31 of the brightest DOGs 
(F 24 /jm >0.8 mjy) in the Bootes Field and with HST imag¬ 
ing. Their sample is dominated by sources with a power law 
in the mid-IR IRAC bands, which is a signature of the pres¬ 
ence of an AGN ( Donley et al.|20'07 |. 

DOGs are believed to be an intermediate AGN phase 
between high-redshift sub millimeter galaxies (SMGs) and 
quasars at z~2 ( Bussmann et al^|2012 l. Accounting for un¬ 
certainty in K which could be as much as a factor of three, 
the median dust mass of our sample is not different from 
those e stimated for hig h-redshift SMGs b y [Magnelli et al.| 
(2012b I (Md ~ 10® Mq). Pope et al. (20081 found that 30% 


of the SMGs from their sample also satisfy the DOG criteria, 
and of those SMG-DOGs, 30 % are AGN dominated. DOGs 
could then be the descendants of these SMGs with similar 
dust content, but representing a more advanced AGN-phase 
than could later quench the star-formation and lead to el¬ 
liptical galaxies. 


6 SUMMARY AND CONCLUSIONS 


We carry out a study that aims to understand the compos¬ 
ite nature of 24/im-bright Dust-Obscured Galaxies (DOGs). 
These sources are a subset of ULIRGs at high redshift 
(z ~ 2) with F 24 /j,m/Fn > 982 . ULIRGs are considered 
to represent an important phase in the evolution of galax¬ 
ies as they are linked to the formation of massive galax¬ 
ies via gas-rich starbursting mergers followed by an AGN- 
driven quenching of the star-formation (e. g., [Sanders et al. 
1988a|b I. Recent studies ( Dey et aL]|2008 Bussmann et al. 
2009 [2012 1 have suggested a similar evolutionary sequence 


where DOGs are an important intermediate phase between 
gas-rich major mergers (traced by submillimeter galaxies, 
SMGs) and quasars at z~2. These studies describe an evo¬ 
lutionary scenario in which the starbursting nature of SMGs 
evolves into the composite nature of DOGs as an underlying 
AGN grows; this is followed by a quasar phase that termi¬ 
nates star formation, leading to the formation of a passive, 
massive elliptical galaxy. Within this context, DOGs could 
provide a key insight to an extremely dusty stage in the 
evolution of galaxies at « ~ 2, where both AGN and star 
formation activity coexist. Their composite nature was un¬ 


til relatively recently inaccessible prior to the availability of 
sensitive mid- to far-infrared data. 

We base our work on a sample of 95 Herschel-detected 
DOG sources. We perform SED-fitting on our sources using 
composite spectra to obtain AGN contributions, dust tem¬ 
peratures and dust masses. We summarize below our results 
and our conclusions: 


(i) DOGs with the brightest 24^m fluxes (F 24 ^m > 
ImJy) present significantly bluer PACS/24/rm colors than 
other 24^m-selected sources. These bluer colors may be ex¬ 
plained by templates containing an AGN contribution of at 
least 25%. 

(ii) Among our sample of 95 sources, 74% are fit by a host 
galaxy template while for 16% require an additional AGN 
component. The remaining 10% of the sample could not be 
properly fit, likely due to inaccurate photometric redshifts. 

(iii) Faint DOG sources with < lO'^® Lq are domi¬ 

nated by star-formation at all redshifts, while DOGs brighter 
than Ls^m > 2 x 10 ^® Lq display a high contribution (> 20 %) 
from an AGN component. 

(iv) DOGs with no significant AGN contribution are 
mainly located within the star-forming main sequence as 
defined in Elbaz et al. (2011 1 . Those identified as AGN- 
DOGs present the lowest IRS (=L/ij/L 8 ) ratio of our sample 
and 50% of them lie below this sequence, with significantly 
lower specific star-formation rates. This results support the 
evolutionary scenario where DOGs may represent a transi¬ 
tion phase between high-redshift starburst-dominated SMGs 
and red-dead ellipticals, passing through an AGN-phase that 
would quench star formation. 

(v) The dust temperature of DOGs peaks at (40 ± 9) K 
and our range of temperatures (24< Td <65 K) is overall in 
good agreement with the liter ature ([Bussmann et al.||2009| 
2012 Melbourne et al. 2012[ Calanog et al. [2013 1. DOGs 


with a contribution from an AGN in the far-IR of at least 
60% have dust temperatures >50 K, suggesting that the 
AGN heats the dust of its host galaxy. We find a median 
dust mass of ~ (3 ±3) x 10® Mq for our sample consistent 
previous analysis in the literature ( [Bussmann et ah ][m^ . 

This work sheds light on DOG sources and their un¬ 
derlying composite nature, bringing unequivocally to light 
that mid-IR bright DOGs are powered by an AGN. The 
submillimeters facilities in the near future, such as GGAT 
and ALMA will provide critical insight to study the AGN 
properties of these obscured ULIRGs at z~2. 

This paper is the first of a series on the panchromatic 
view of DOGs. In this paper we focus on their far-IR prop¬ 
erties, while in the upcoming papers we will focus on their 
X-ray properties based on X-ray stacking analysis and on 
their contribution to the Cosmic X-ray Background. The 
far-IR/radio correlation of these sources and their radio 
properties using JVLA-COSMOS observations presented in 


Smolcic et al. (2014 \ will be detailed in a forthcoming paper. 
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Table 3. DOGs sample 


DOG ID 

redshift 

F 24 

(mJy) 

Fioo 

(mJy) 

F 16 O 

(mJy) 

F 250 

(mJy) 

F 350 

(mJy) 

F 500 

(mJy) 

flag AGN 


Lir 

(L©) 

0 

1.87 

0.386 

9.798 

14.33 

35.05 

17.40 

14.67 

0 

2.90 

1.02e-|-12 

1 

1.60 

0.178 

11.48 

10.53 

23.54 

16.14 

10.17 

0 

9.59 

8.41e-|-ll 

2 

2.65 

0.176 

9.251 

11.81 

21.27 

19.14 

14.26 

0 

6.22 

3.41e-|-12 

3 

2.34 

0.491 

16.73 

26.67 

27.36 

20.72 

7.077 

0 

3.00 

2.48e-|-12 

4 

3.00 

0.097 

9.961 

27.31 

25.93 

14.49 

3.826 

3 

0.409 

3.67e-|-12 

5 

1.00 

0.404 

9.088 

20.96 

12.49 

8.223 

6.096 

0 

37.8 

1.30e-|-12 

6 

1.40 

0.440 

7.806 

16.88 

15.89 

13.25 

9.133 

0 

20.5 

6.84e-|-ll 

7 

2.85 

0.326 

8.145 

13.69 

22.02 

18.26 

5.261 

2 

2.19 

1.20e-|-12 

8 

2.34 

0.402 

8.467 

24.79 

22.32 

20.23 

9.081 

0 

3.00 

2.48e-|-12 

9 

1.89 

0.312 

15.11 

10.92 

10.72 

10.60 

8.063 

0 

7.78 

1.49e-|-12 

10 

1.14 

0.103 

5.576 

16.27 

15.17 

17.43 

5.545 

0 

3.56 

4.15e-|-ll 

11 

1.98 

0.156 

14.28 

18.67 

25.49 

17.19 

2.676 

-99 

-99 

-99 

12 

1.88 

0.236 

11.65 

29.26 

19.52 

7.441 

3.868 

0 

24.5 

1.587e-|-12 

13 

1.81 

0.444 

14.62 

43.96 

31.01 

17.75 

7.619 

0 

0.0739 

9.502e-|-ll 

14 

1.79 

0.201 

6.918 

32.55 

30.33 

19.84 

9.915 

0 

31.7 

2.460e-|-12 

15 

2.30 

0.326 

8.989 

18.93 

30.08 

21.45 

20.49 

0 

8.01 

2.068e-|-12 

16 

2.04 

0.656 

9.732 

20.43 

30.88 

34.16 

15.04 

4 

0.668 

1.753e-|-12 

17 

2.34 

0.716 

11.18 

20.34 

20.45 

16.57 

8.984 

0 

0.311 

1.760e-|-12 

18 

2.88 

0.330 

17.99 

54.18 

43.51 

35.79 

32.05 

0 

12.6 

4.012e-|-12 

19 

2.74 

0.223 

12.73 

22.93 

32.10 

20.40 

5.487 

0 

0.509 

4.354e-|-12 

20 

2.41 

0.665 

16.44 

27.88 

17.78 

14.48 

0.769 

0 

2.57 

1.842e-|-12 

21 

1.90 

0.363 

6.506 

9.178 

21.98 

19.80 

13.77 

0 

4.31 

1.538e-|-12 

22 

2.15 

0.319 

9.573 

31.35 

50.37 

43.09 

29.44 

0 

29.9 

1.798e-|-12 

23 

2.11 

0.400 

9.302 

26.70 

34.57 

26.41 

13.75 

0 

15.8 

2.122e-|-12 

24 

2.14 

0.208 

5.767 

15.81 

18.79 

12.01 

5.824 

0 

32.7 

2.412e-|-12 

25 

2.50 

0.414 

4.650 

14.08 

35.78 

28.10 

10.19 

-99 

-99 

-99 

26 

2.29 

0.485 

7.434 

15.55 

26.60 

25.71 

25.31 

0 

9.11 

3.31e-|-12 

27 

1.84 

0.285 

7.054 

23.17 

24.96 

21.93 

9.473 

0 

5.71 

9.19e-|-ll 

28 

1.79 

0.400 

7.824 

16.42 

12.27 

14.91 

6.783 

0 

3.61 

1.59e-|-12 

29 

1.80 

0.532 

7.399 

21.88 

28.97 

24.07 

13.26 

-99 

-99 

-99 

30 

1.91 

0.287 

8.802 

23.61 

27.99 

26.50 

17.52 

0 

2.37 

1.14e-|-12 

31 

1.80 

0.601 

10.28 

32.24 

50.22 

57.26 

37.58 

0 

0.503 

1.67e-|-12 

32 

1.43 

0.271 

27.75 

50.52 

54.32 

22.67 

22.84 

0 

47.8 

1.42e-|-12 

33 

1.24 

0.405 

7.274 

12.76 

24.40 

25.07 

12.53 

0 

0.0981 

6.33e-|-ll 

34 

1.73 

0.090 

3.589 

13.36 

18.04 

8.671 

1.258 

0 

1.04 

1.04e-|-12 

35 

1.83 

0.669 

7.372 

18.38 

29.18 

21.97 

15.27 

0 

3.06 

1.73e-|-12 

36 

1.88 

0.388 

12.93 

17.25 

11.61 

5.563 

3.736 

0 

13.3 

1.63e-|-12 

37 

2.18 

0.185 

11.57 

16.36 

20.27 

11.07 

7.944 

0 

4.96 

2.85e-|-12 

38 

1.68 

0.846 

7.959 

10.58 

43.39 

42.47 

40.47 

0 

2.16 

9.40e-|-ll 

39 

1.56 

0.280 

9.949 

15.29 

16.03 

18.52 

12.97 

0 

10.0 

1.25e-|-12 

40 

1.71 

0.412 

17.15 

40.69 

38.81 

23.32 

9.754 

0 

4.66 

1.02e-|-12 

41 

2.00 

0.737 

10.41 

21.83 

31.49 

28.50 

17.59 

0 

9.79 

1.17e-|-12 

42 

1.79 

0.282 

5.335 

14.58 

22.06 

15.15 

3.131 

0 

8.24 

1.37e-|-12 

43 

2.75 

0.277 

7.221 

21.95 

14.15 

16.22 

13.70 

3 

15.7 

8.33e-|-12 

44 

2.03 

1.038 

29.37 

43.28 

34.96 

28.66 

7.363 

2 

13.4 

4.74e-|-ll 

45 

1.91 

0.248 

5.540 

14.59 

26.78 

19.12 

7.714 

0 

5.65 

1.43e-|-12 

46 

1.81 

0.551 

16.28 

28.14 

17.47 

10.42 

6.205 

0 

4.42 

1.72e-|-12 

47 

2.12 

0.179 

5.090 

8.722 

25.19 

24.05 

19.96 

-99 

-99 

-99 

48 

2.35 

0.386 

9.840 

30.80 

33.42 

44.45 

28.70 

0 

11.8 

4.28e-|-12 

49 

1.96 

0.130 

6.263 

12.26 

20.61 

23.57 

12.48 

0 

1.40 

1.48e-|-12 

50 

2.92 

0.296 

9.159 

25.80 

40.72 

28.25 

21.91 

0 

8.62 

3.72e-|-12 

51 

2.55 

0.559 

9.291 

20.90 

20.95 

4.990 

16.31 

0 

1.26 

1.82e-|-12 

52 

1.94 

0.208 

6.285 

11.61 

19.85 

18.14 

7.037 

0 

2.37 

1.72e-|-12 

53 

2.26 

0.379 

9.616 

28.90 

38.44 

34.29 

20.98 

0 

8.56 

1.51e-|-12 

54 

1.88 

0.704 

4.882 

29.71 

30.02 

36.79 

7.011 

0 

0.672 

1.41e-|-12 

55 

1.61 

0.187 

7.526 

18.13 

19.77 

23.50 

15.58 

0 

25.1 

2.08e-|-12 

56 

1.61 

0.284 

7.408 

19.38 

21.56 

17.52 

9.390 

0 

2.67 

l.OOe-l-12 

57 

1.73 

0.537 

5.251 

16.61 

29.93 

25.76 

9.697 

0 

24.1 

1.73e-|-12 

58 

2.77 

0.165 

7.150 

14.23 

13.25 

8.367 

8.701 

2 

20.9 

2.33e-|-12 

59 

1.94 

0.330 

14.65 

24.56 

26.93 

21.13 

10.50 

0 

0.369 

1.94e-|-12 

60 

1.27 

0.433 

10.52 

27.31 

44.89 

44.77 

34.52 

0 

9.25 

7.57e-|-ll 


flag AGN is the contribution from an AGN to the host galaxy obtained from DecompIR: (0): only host galaxy, (1): % AGN ^ 20%, (2): 

20 < % AGN ^ 40%, (3): 40 < % AGN ^ 70%, (4): % AGN > 70% 
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Table 3 — continued 


DOG ID 

redshift 

F 24 

(mjy) 

Fioo 

(mJy) 

F 16 O 

(mJy) 

F 250 

(mJy) 

F 350 

(mJy) 

F 500 

(mJy) 

flag AGN 


(L©) 

61 

2.00 

0.745 

15.31 

25.86 

38.35 

47.96 

29.10 

0 

50.1 

3.18e-|-12 

62 

1.58 

0.505 

6.231 

31.96 

8.530 

45.61 

8.272 

0 

7.96 

1.22e-|-12 

63 

1.88 

0.621 

17.10 

39.68 

41.01 

29.16 

21.66 

0 

6.60 

6.77e-|-ll 

64 

2.55 

0.366 

8.821 

15.62 

8.706 

3.102 

6.735 

0 

2.24 

2.61e-|-12 

65 

1.76 

0.419 

9.343 

19.55 

28.59 

27.63 

20.37 

-99 

-99 

-99 

66 

2.03 

0.407 

14.68 

28.11 

22.86 

15.96 

10.10 

0 

2.44 

2.69e-|-12 

67 

2.85 

0.497 

16.95 

27.73 

10.42 

8.087 

11.99 

0 

14.7 

4.00e-|-12 

68 

1.94 

0.943 

13.47 

26.86 

15.43 

7.120 

2.860 

-99 

-99 

-99 

69 

1.62 

0.925 

43.88 

66.54 

54.21 

32.34 

6.816 

-99 

-99 

-99 

70 

1.93 

0.437 

7.177 

19.51 

21.00 

8.710 

8.849 

0 

1.69 

1.39e-|-12 

71 

2.53 

0.132 

16.52 

44.25 

51.89 

34.34 

23.83 

3 

1.20 

2.98e-|-12 

72 

2.91 

0.555 

7.702 

23.22 

19.01 

19.22 

10.53 

-99 

-99 

-99 

73 

1.61 

0.395 

13.92 

23.40 

31.45 

13.30 

11.82 

4 

1.44 

2.49e-|-12 

74 

1.61 

0.266 

7.749 

20.32 

20.94 

10.23 

1.243 

3 

0.200 

3.54e-|-ll 

75 

2.70 

0.256 

6.372 

13.61 

29.14 

24.31 

23.73 

0 

1.67 

3.36e-|-12 

76 

2.92 

0.465 

13.72 

41.37 

39.27 

37.63 

28.69 

0 

1.60 

2.94e-|-12 

77 

2.00 

1.487 

10.00 

17.00 

29.15 

27.47 

5.519 

0 

14.1 

1.48e-|-12 

78 

1.98 

0.554 

22.49 

71.89 

74.66 

50.33 

50.14 

-99 

-99 

-99 

79 

2.11 

0.359 

14.87 

21.01 

25.52 

18.83 

15.39 

0 

11.0 

2.55e-|-12 

80 

1.89 

4.742 

20.84 

26.75 

20.47 

16.42 

8.412 

2 

0.761 

6.38e-|-ll 

81 

2.33 

1.385 

30.01 

59.63 

62.33 

50.04 

24.83 

0 

3.97 

1.31e-|-12 

82 

1.97 

0.625 

12.46 

20.73 

26.62 

27.96 

11.84 

0 

5.23 

1.26e-|-12 

83 

1.98 

1.059 

28.67 

36.87 

32.63 

40.06 

27.67 

3 

15.6 

2.07e-|-12 

84 

1.91 

0.578 

7.835 

13.90 

48.09 

51.41 

32.16 

0 

3.45 

1.74e-|-12 

85 

2.58 

1.917 

14.49 

18.95 

13.09 

9.595 

2.781 

0 

15.9 

2.26e-|-12 

86 

2.91 

0.759 

11.60 

21.18 

35.11 

49.91 

24.43 

2 

0.021 

l.OOe-l-12 

87 

2.34 

0.637 

7.290 

17.49 

20.34 

23.81 

13.98 

0 

6.19 

2.14e-|-12 

88 

2.64 

3.744 

55.74 

102.9 

100.3 

59.92 

55.24 

1 

8.16 

3.29e-|-12 

89 

2.20 

0.860 

12.87 

22.32 

12.75 

19.16 

21.22 

0 

19.6 

5.25e-|-12 

90 

2.89 

0.671 

16.75 

27.59 

41.20 

36.35 

4.689 

1 

4.33 

5.54e-|-ll 

91 

2.88 

0.931 

5.011 

11.63 

15.30 

13.39 

16.10 

1 

5.19 

8.52e-|-ll 

92 

2.60 

2.392 

17.72 

34.43 

49.46 

32.61 

10.89 

0 

8.08 

2.23e-|-12 

93 

1.73 

3.131 

8.783 

10.84 

18.91 

13.82 

2.563 

0 

12.8 

1.76e-|-12 

94 

1.75 

1.559 

20.80 

20.48 

16.66 

8.258 

4.248 

0 

6.89 

2.05e-|-12 
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